Sounding the Solar Cycle with 
Helioseismology: Implications for 
Asteroseismology 



By W. J. Chaplin 



1. Introduction 

My brief for the IAC Winter School was to cover observational results on helioseismol- 
ogy, flagging where possible implications of those results for the asteroseismic study of 
solar-type stars. My desire to make such links meant that I concentrated largely upon 
results for low angular-degree (low-/) solar p modes, in particular results derived from 
"Sun-as-a-star" observations (which are of course most instructive for the transfer of ex- 
perience from helioseismology to asteroseismology). The lectures covered many aspects 
of helioseismology - modern helioseismology is a diverse field. In these notes, rather than 
discuss each aspect to a moderate level of detail, I have instead made the decision to 
concentrate upon one theme, that of "sounding" the solar activity cycle with helioseis- 
mology. I cover the topics from the lectures and I also include some new material, relating 
both to the lecture topics and other aspects I did not have time to cover. Implications 
for asteroseismology are developed and discussed throughout. 

The availability of long timeseries data on solar-type stars, courtesy of the NASA 
Kepler Mission (Gilliland et al. 2010; Chaplin et al. 2010) and the French-led CoRoT 
satellite (Appourchaux et al. 2008), is now making it possible to "sound" stellar cycles 
with asteroseismology. The prospects for such studies have been considered in some depth 
(e.g., Chaplin et al. 2007, 2008a; Metcalfe et al. 2008; Karoff et al. 2009), and in the last 
year the first convincing results on stellar-cycle variations of the p-mode frequencies 
of a solar- type star (the F-type star HD49933) were reported by Garcia et al. (2010). 
This result is important for two reasons: first, the obvious one of being the first such 
result, thereby demonstrating the feasibility of such studies; and second, the period of 
the stellar cycle was evidently significantly shorter than the 11-yr period of the Sun 
(probably between 1 and 2yr). If other similar stars show similar short- length cycles, 
there is the prospect of being able to "sound" perhaps two or more complete cycles of 
such stars with Kepler (assuming the mission is extended, as expected, to 6.5 yr or more). 
The results on HD49933 may be consistent with the paradigm that stars divide into two 
groups, activity-wise, with stars in each group displaying a similar number of rotation 
periods per cycle period (e.g., see Bohm-Vitense 2007), meaning solar- type stars with 
short rotation periods - HD49933 has a surface rotation period of about 3 days - tend to 
have short cycle periods. We note that Metcalfe et al. (2010) recently found another F- 
type star with a short (1.6 yr) cycle period (using chromospheric H & K data). Extension 
of the Kepler Mission will of course also open the possibility of detecting full swings in 
activity in stars with cycles having periods up to approximately the length of the solar 
cycle. 

The rest of my notes break down as follows. Section [5] gives an introductory overview 
of the solar cycle, as seen in helioseismic data. A brief history of observations of solar 
cycle changes in low-angular degree (low-Z) solar p modes is given in Section 12.11 Then, 
in Section [2~2l we consider the causes of the observed changes in the mode frequencies; 
and in Section 12. 3[ we discuss variations in the mode powers and damping rates, and 
what the relative sizes of those changes imply for the underlying cause. 
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Section [3] considers several subtle ways in which the stellar activity cycles can affect 
values of, and inferences made from, the mode parameters. Concepts are introduced using 
the example of the solar cycle, and the impact it has on low-/ p modes observed in Sun- 
as-a-star data. Implications for asteroseismic observations of solar-type stars are then 
developed. We start in Section [3711 with a discussion of the impact of stellar activity on 
estimates of the mode frequencies. This is followed in Section |3~21 by a similar discussion 
for frequency separation ratios. Finally, Section 13.31 explains how mode peaks in the 
frequency-power spectrum can be "distorted" by stellar cycles, rendering commonly used 
fitting models inappropriate. 

In Scction|4]we develop a simple model to illustrate the impact on the mode frequencies 
of the range in latitudes covered by near-surface magnetic activity on solar-type stars, 
and discuss how the angle of inclination can affect significantly the observed frequency 
shifts (because that angle affects which mode components are visible in the observations). 
We also show how measurements of the frequency shifts of modes having different angular 
and azimuthal degrees may be used to make inference on the spatial distribution of the 
near-surface magnetic activity on solar-type stars. 

We end in Section [5] by thinking somewhat longer-term, and consider how much low-/ 
data would be needed to measure evolutionary changes of the solar p-mode frequencies. 

2. The seismic solar cycle: overview 

A rich, and diverse, body of observational data is now available on temporal varia- 
tions of the properties of the global solar p modes. The signatures of these variations are 
correlated strongly with the well-known 11-year cycle of surface activity. The search for 
temporal variations of the p-mode properties began in the early 1980s, following accu- 
mulation of several years of global seismic data. The first positive result was reported 
by Woodard & Noyes (1985), who found evidence in observations made by the Active 
Cavity Radiometer Irradiance Monitor (ACRIM) instrument, on board the Solar Max- 
imum Mission (SMM) satellite, for a systematic decrease of the frequencies of low-/ p 
modes between 1980 and 1984. The first year coincided with high levels of global sur- 
face activity, while during the latter period activity levels were much lower. The modes 
appeared to be responding to the Sun's 11-year cycle of magnetic activity. Woodard & 
Noyes found that the frequencies of the most prominent modes had decreased by roughly 
1 part in 10 000 between the activity maximum and minimum of the cycle. By the late 
1980s, an in-depth study of frequency variations of global p modes, observed in the Big 
Bear data, had demonstrated that the agent of change was confined to the outer layers 
of the interior (Libbrecht & Woodard, 1990). 

Accumulation of data from the new networks and instruments has made it possible 
to study the frequency variations to unprecedented levels of detail, and has revealed 
signatures of subtle, structural change in the sub-surface layers. The discovery of solar- 
cycle variations in mode parameters associated with the excitation and damping (e.g., 
power, damping rate, and peak asymmetry) followed. Patterns of flow that penetrate a 
substantial fraction of the convection zone have also been uncovered as well as possibly 
signatures of changes in the rotation rate of the layers that straddle the tachocline, 
and much more recently evidence for a quasi-periodic 2-yr signal, superimposed on the 
solar-cycle variations of the mode frequencies. 

The modern seismic data give unprecedented precision on measurements of frequency 
shifts. Examples of average frequency shifts for low-/ data are shown in Fig. Q] for 
Sun-as-a-star data collected by the ground-based Birmingham Solar-Oscillations Net- 
work (BiSON), and Global Oscillations at Low Frequency (GOLF) instrument on board 
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Figure 1. Average frequency shifts of most prominent low-/ solar p modes, as measured in 
BiSON data (top panel) and GOLF data (bottom panel) over the last two 11-yr activity cycles. 
The three curves show results for averages made over different ranges in frequency: 1880 to 
3710 /iHz (diamonds, joined by solid line); 1880 to 2770 /iHz (crosses, joined by dotted line); 
2820 to 3710 /xHz (triangles, joined by dashed line). (From Fletcher et al. 2010.) 



the ESA/NASA Solar and Heliospheric Observatory (SOHO). From observations of the 
medium-Z frequency shifts it is possible to produce surface maps showing the strength of 
the solar-cycle shifts as a function of latitude and time (Howe, Komm, & Hill, 2002), like 
the example shown in Figure [5] (which is made from Global Oscillations Network Group 
(GONG) data). These maps bear a striking resemblance to the butterfly diagrams that 
show spatial variations in the strength of the surface magnetic field over time. The im- 
plication is that the frequency shift of a given mode depends on the strength of that 
component of the surface magnetic field that has the same spherical harmonic projection 
on the surface. Similar maps may also be made for variations observed in the p-mode 
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Figure 2. Mode frequency shifts (in iiHz) as a function of time and latitude. The values come 
from analysis of GONG data. The contour lines indicate the surface magnetic activity. (Figure 
courtesy of R. Howe.) 

powers and damping rates (Komm, Howe, & Hill, 2002), which, like the frequency maps, 
show a close spatial and temporal correspondence with the evolution of active-region field 
(Figure©. 

2.1. The seismic solar cycle at low angular degree 

As noted above, the first evidence for activity-related changes to the low-? p modes of the 
Sun was reported by Woodard & Noyes (1985). These changes were soon confirmed (and 
the results extended) by Palle ct al. (1989) and Elsworth et al. (1990). Almost two decades 
on, there now exists an extensive literature devoted to studies of low-Z mode parameter 
variations. Since the frequencies are by far the most precisely determined parameters, it 
is hardly surprising that analyses of the frequencies threw up the first positive results 
for solar-cycle changes. Evidence for changes in mode power followed next (Palle et al. 
1990a; Anguera-Gubau et al. 1992; Elsworth et al. 1994). The first tentative claims for 
changes in mode linewidth (damping) were made by Palle et al. (1990b). Subsequently, 
Toutain & Wehrli (1997) and Appourchaux (1998) provided stronger evidence in support 
of an increase in damping with activity, and these claims were confirmed beyond all 
doubt by Chaplin et al. (2000). (Komm, Howe & Hill (2000) did likewise for medium 
I modes at about the same time.) Peak asymmetry is the most recent addition to the 
list of parameters that show solar-cycle variations (Jimenez-Reyes et al. 2007). Careful 
measurement of variations in the powers, damping rates and peak asymmetries - all 
parameters associated with the excitation and damping - are allowing studies to be 
made of the impact of the solar cycle on the convection properties in the near-surface 
layers. 

Frequencies of modes below w 4000 /iHz are observed to increase with increasing ac- 



W. J. Chaplin: Sounding the Solar Cycle 



5 




1996 1998 2000 2002 2004 

Date (years) 



Figure 3. Mode linewidth (in /iHz) as a function of time and latitude. The values come from 
analysis of GONG data. The contour lines indicate the surface magnetic activity. (Courtesy of 
R. Howe.) 



tivity. For the most prominent modes at around 3000 fiHz, the size of the shift is about 
0.4 /iHz between activity minimum and maximum. Furthermore, higher frequency modes 
experience a larger shift than their lower-frequency counterparts. This frequency depen- 
dence suggests that the perturbations responsible for the frequency shifts are located 
very close to the solar surface. The upper turning points of the modes - which for low- 
degrec modes arc effectively independent of I - lie deeper in the Sun for low-frequency 
modes than they do for high-frequency modes: higher-frequency modes are as such more 
sensitive to surface perturbations. 

That the shifts do not scale like v n ijL, where v n i is mode frequency and L = sjl{l + 1), 
rules out the possibility that the perturbation is spread throughout a significant fraction 
of the solar interior. This may be understood by thinking classically, in terms of ray 
paths followed by the acoustic waves. When a wave reaches the lower turning point of its 
cavity it will by definition be moving horizontally, and its phase speed will be equal to 

u n i 2ixv rd R 

c = ~ = ~JW+T) a (2J) 

where k is the horizontal wavenumber and R the outer cavity radius. The ratio v n i/L 
therefore maps to the location of the lower boundary of the cavity, and hence the cavity 
size. Since the shifts do not scale like this ratio, we conclude that the perturbation must 
be confined within a narrow layer, and not spread so widely that it covers the entirety 
of the cavities of many of the modes. 

At frequencies above 4000 ^Hz the size of the shift decreases, and also changes sign 
above ~ 4500 ^Hz meaning that these very high-frequency modes suffer a reduction in 
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frequency as activity levels rise (Anguera-Gubau et al. 1992; Chaplin et al. 1998; Jimenez- 
Reyes et al. 2001; Gelly et al. 2002; Salabert et al. 2004). 

Detailed comparison of the \ow-l frequency shifts with changes in various disc-averaged 
proxies of global surface activity provides further tangible input to the solar cycle studies. 
This is because different proxies show differing sensitivity to various components of the 
surface activity. While the changes in frequency are observed to correlate fairly well with 
contemporaneous changes in global proxies, the match is far from perfect. Jimenez- Reyes 
et al. (1998) were the first to show that the relation of the frequency shifts to variations 
in the proxies was markedly different on the rising and falling parts of the 11-yr Schwabe 
cycle. Since the magnitudes of the shifts should reflect the different spatial sensitivities of 
modes of different angular and azimuthal degree to the time-dependent variation of the 
surface distribution of the activity, a better choice for the activity proxy would clearly be 
one that has been decomposed to have a similar spatial distribution as the mode under 
study. Chaplin et al. (2004a) and Jimenez-Reyes et al. (2004) have shown that the sizes of 
the \ow-l shifts do indeed scale better with activity proxies that have the same spherical 
harmonic projection as the modes. 

Chaplin et al., (2007b) compared frequency changes in 30 years of BiSON data with 
variations in six well-known activity proxies. Interestingly, they found that only activity 
proxies having good sensitivity to the effects of weak-component magnetic flux - which 
is more widely distributed in latitude than the strong flux in the active regions - were 
able to follow the frequency shifts consistently over the three cycles. 

The unusual behaviour during the most recent solar minimum (straddling solar cycles 
23 and 24) of many diagnostics and probes of solar activity has raised considerable 
interest and debate in the scientific community (e.g., see the summary by Sheeley 2010). 
The minimum was unusually, and unexpectedly, extended and deep. Polar magnetic fields 
were very weak, and the open flux was diminished compared to other preceding minima. 

Helioseismology has been used to probe the behaviour of sub-surface flows during 
the minimum. Howe et al. (2009) found that the equatorward progression of the lower 
branches of the so-called torsional oscillations (east-west flows) was late in starting com- 
pared to previous cycles. They flagged this delayed migration as a possible pre-cursor of 
the delayed onset of cycle 24. The meridional (north-south) flow also carries a signature 
of the solar cycle, which converges toward the active-region latitudes and also intensifies 
in strength as activity increases: Gonzalez-Hernandez (2010) found that during the cur- 
rent minimum this component had developed to detectable levels even before the visual 
onset of magnetic activity on the solar surface. 

The globally coherent acoustic properties of the recent solar minimum have been stud- 
ied extensively with low-degree p modes (by Broomhall et al. 2009 and Salabert et al. 
2009) and medium-degree p modes (Tripathy et al. 2010). These studies have shown that 
while the surface proxies of activity (e.g., the 10.7-cm radio flux) were quiescent and 
very stable during the minimum, the p-mode frequencies showed much more variability. 
Tripathy et al. (2010) noted further surprising behaviour compared to the previous cy- 
cle 22-23 minimum, i.e., an apparent anti-correlation of the p-mode frequency shifts and 
the surface proxies of activity. 

Broomhall et al. (2009) had suggested the possible presence of a quasi-biennial modu- 
lation of the frequencies of the low-degree modes, superimposed upon the well-established 
~ 11-yr variation of the frequencies. This has since been confirmed by further in-depth 
analysis, which reveals a signature that is consistent in the frequencies extracted from Bi- 
SON and GOLF data (Fletcher et al. 2010), as shown in Fig.QJ This plots the frequency 
residuals that remain after the long-term solar-cycle variation has been removed from 
the average frequency shifts. The residuals show significant variability, with a period of 
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Figure 4. Frequency residuals that remain after the long-term solar-cycle variation has been 
removed from the average frequency shifts of BiSON data (top panel) and GOLF data (bottom 
panel). The three curves in each panel show results for averages made over different ranges in 
frequency: 1880 to 3710 /^Hz (diamonds, joined by solid line); 1880 to 2770 /iHz (crosses, joined 
by dotted line); 2820 to 3710 /iHz (triangles, joined by dashed line). (Plot from Fletcher et al. 
2009.) 



around 2 yr, that variability being more pronounced at times of higher surface activity. 
The fact that this biennial signature has similar amplitude in the low-frequency and the 
high-frequency modes used in this analysis suggests that its origins lie deeper than the 
very superficial layers responsible for the 11-year shifts. 

Finally in this section on the frequencies, we note that from appropriate combinations 
of low-i frequencies Verner et al. (2006) uncovered apparent solar-cycle variations in the 
amplitude of the depression in the adiabatic index, T\, in the He II zone. These variations 
presumably reflect the impact of the changing activity on the equation of state of the 
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gas in the layer, and confirmed the findings of Basu and Mandel (2004), who used the 
more numerous data available from medium-? frequencies. 

2.2. What is the cause of the frequency shifts? 

Broadly speaking, the magnetic fields can affect the modes in two ways. They can do 
so directly, by the action of the Lorentz force on the gas. This provides an additional 
restoring force, the result being an increase of frequency, and the appearance of new 
modes. Magnetic fields can also influence matters indirectly, by affecting the physical 
properties in the mode cavities and, as a result, the propagation of the acoustic waves 
within them. This indirect effect can act both ways, to either increase or decrease the 
frequencies. 

We begin this section with a back-of-the-envelope calculation, which makes (an admit- 
tedly) approximate prediction of the impact on the frequency of a typical low-? mode of 
changes in stratification brought about by the near-surface magnetic field in sunspots. 
The calculation is instructive, in that the size of the predicted frequency shift it gives is 
broadly in line with the observations, providing further evidence in support of the shifts 
being caused by near-surface perturbations due to the presence of magnetic field. 

As we have seen, frequencies of the most prominent modes increase with increasing 
magnetic activity. Here, we consider the indirect effect of the near-surface magnetic field 
on the near-surface properties, and hence the frequencies of these modes. If magnetic 
fields modify the surface properties in such a way as to reduce the effective size of the 
mode cavity, the required increase of frequency will result. 

In regions of strong magnetic field - such as those occupied by sunspots - there will 
be a gas pressure deficit (assuming those regions to be in pressure equilibrium with their 
field- free surroundings). This is because gas and magnetic pressure combine within the 
magnetic regions, while only gas pressure acts in the field-free regions. Sunspots are also 
characterized by a reduced temperature, relative to the surroundings. The central part of 
the spot therefore exhibits a lower pressure, temperature and density than the surround- 
ings, resulting in the so-called Wilson Depression. Values of pressure, temperature and 
density found at the surface in the field-free plasma are only reached at some depth be- 
neath the surface in strong-field regions. Recent measurements suggest that for a sunspot 
the typical size of this depression is about 1000 km (Watson et al. 2009). 

Let us assume that 1000 km corresponds approximately to the amount, 5R, by which 
the mode cavities are reduced in size beneath sunspots. The fractional area of the solar 
surface occupied by sunspots reaches ~ 0.5% at modern cycle maxima. We therefore 
obtain a net, surface-averaged estimate of SR via 

(SR) w 0.5 % x 1000 km « 5 km. 

The sound speed, c, at the solar surface is approximately lOkrns -1 , implying a reduction 
in the travel time in the mode cavity, due to the "shrinkage" at the surface, of 

STm Skm/lOkms" 1 w 0.5s. 

The travel time across the cavity, T, is related to the large frequency separation, Au, 
via: 

Ai/~^2jf*^ ~ (2T)- 1 . (2.2) 

For the Sun, Au = 135 ^Hz, implying T rts 3700 s. The fractional change (reduction) in 
T is therefore: 

ST/T ~ 0.5/3700 ~ 1.4 x 10~ 4 . 
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Since 5v/v = —6T/T, the predicted increase in frequency of a mode at s» 3000 /iHz will 
be 

5v = 3000 x 1.4 x 10~ 4 ~ 0.4/zHz. 

This estimate is of a very similar size to the observed frequency shifts. 

What does detailed modelling suggest? Perhaps the most significant contribution of 
recent years in this area is that of Dziembowski and Goode (2005). Their results suggest 
that the indirect effects dominate the perturbations, and that the magnetic fields are 
too weak in the near-surface layers for the direct effect to contribute significantly to the 
observed frequency shifts. However, Dziembowski and Goode also found some evidence to 
suggest that the direct effect may play a more important role for low-frequency modes, at 
depths beneath the surface where the magnetic field is strong enough to give a significant 
direct contribution to the frequency shifts. We now go on to explain how dependence of 
the frequency shifts on mode inertia and mode frequency can tell us something about 
the location and nature of the perturbations. 

We begin by noting that when the frequency shifts are multiplied by the mode inertia, 
and then normalized by the inertia of a radial mode of the same frequency, the mod- 
ified shifts are found to be a function of frequency alone. This in effect removes any I 
dependence of the shifts (at fixed frequency, the higher the Z, the larger is the observed 
frequency shift). The observed I dependence may be understood in terms of, for example, 
a physical interpretation of the mode inertia. The normalized inertia, I n i , may be defined 
according to (Christensen-Dalsgaard & Berthomieu 1991): 

Im = Mq 1 f = 47TM- 1 / \Z\ 2 pr 2 dr = M nl /M Q (2.3) 

Jv Jo 

where £ is the (surface-normalized) displacement associated with the mode, and the 
integration is performed over the volume V of the Sun, which has mass Mq. The mode 
mass M n i is therefore the interior mass affected by the perturbations associated with 
the mode. As Z increases, so M n i decreases, and the more sensitive a mode will be to 
a near-surface perturbation of a given size (giving a larger frequency shift). One may 
therefore render the shifts Z independent by multiplying them by the inertia ratio Q n i 
(Christensen-Dalsgaard & Berthomieu 1991), which is given by: 

Qnl = Inl/I{Vnl) (2.4) 

Multiplication of the raw shifts by Q„i is indeed seen to collapse the shifts of different Z 
onto a single curve. As shown in Fig. [5l this then allows one to combine data spanning 
a range in I, which reduces errors, giving tighter constraints on frequency dependence of 
the shifts. 

Chaplin et al. (2001) studied in detail the frequency dependence of the inertia-ratio- 
corrected shifts, 5iy n iQ n i, of both low-Z modes and medium-Z modes up to I — 150. They 
fitted these data to a power law of the form 

SVnlQnl = -^Ku ( 2 -5) 

where the power-law index is a and c is a constant. They found that a ~ 2 for v > 
2500 /iHz, while a is approximately zero for v < 2500 /xHz. Rabello-Soares et al. (2008) 
repeated the exercise for medium-Z and high-Z modes. They extracted very similar be- 
haviour, and, because they had more medium- and high-Z data at low frequencies, they 
were able to extend the analysis to v < 2000 /xHz where they found that a appears to 
change sign and go negative. The high-Z modes provide important information since they 
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Figure 5. Top panel: Frequency shifts, Sv, plotted as a function of mode frequency, for results 
from different instruments over different ranges in I (see plot annotation). Bottom panel: Scaled 
frequency shifts, after multiplication by the inertia ratio Q ni . (From Chaplin et al. 2001.) 



are confined in the layers close to the surface where the physical changes responsible for 
the frequency shifts are also located. 

What do the values of a imply? If the perturbation is located within the photosphere 
(but is confined to extend over an extent no longer than one pressure scale height), 
then we might expect a ~ 3 (e.g., Libbrecht & Woodard 1990). If instead the pertur- 
bation extends beneath the surface, the frequency dependence will be weaker and a will 
get smaller (Gough 1990). The observed values of a in the above mentioned frequency 
regimes are both less than 3, suggesting that the perturbations cannot arise solely in 
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the photosphere. Since a is smaller for the lower-frequency modes, this suggests that the 
perturbation extends to greater depths the lower in frequency one goes. This is consistent 
with the inferences made by Goode & Dziembowski (2005). 

2.3. Solar-cycle variations of mode power and mode damping 

Variations in mode damping have now been uncovered in all the major low-Z datasets 
(see Chaplin 2004 and references therein). The trend found is an increase of about 20 per 
cent between activity minimum and maximum, with some suggestion of the variations 
being peaked in size at about ~ 3000 /iHz. Mode powers are at the same time observed to 
decrease by about the same fractional amount, while the mode heights decrease by twice 
the amount. Use of the analogy of a damped, randomly forced oscillator is particularly 
instructive for understanding these results. 

We define our oscillator according to the equation: 

x(t) + 2r)x(t) + (2Trv ) 2 x(t) = KS(t - t ), (2.6) 

where x(t) is the displacement, vq the natural frequency of the oscillator, r\ the linear 
damping constant, and K is the amplitude of the forcing function (assumed to be a 
random Gaussian variable), with "kicks" applied at times t n , 5(t — to) being the delta 
function. 

Provided that T{y) - the frequency spectrum of the forcing function - is a slowly- 
varying function of v, and r\ << 2-kvq, the power spectral density (PSD) in the frequency 
domain will be a Lorentzian profile, i.e., 



PSD(i/) oc T(v) 1 + — — ± . (2.7) 




This holds for both the spectrum of the displacement, x(t), and the velocity, x(t). The 
FWHM of the peak in cyclic frequency is 

A=^. (2.8) 

The maximum power spectral density, or height, of the resonant peak is: 

ffoc^M. (2.9) 
r, 2 

The total mean-square power (variance in the time domain) is proportional to peak height 
times width, i.e., P cx HA, so that 

Poc^M. (2.10) 
V 

The energy (kinetic plus potential) of a resonant mode with associated inertia / is given 
by: 

E = MP. (2.11) 

The rate at which energy is supplied to (and dissipated by) the modes, dE/dt, is readily 
derived by again having recourse to the oscillator analogy. The amplitude of the oscillator 
is attenuated in time by the factor exp (— r/i), and its energy is proportional to the 
amplitude squared. Hence, we may write 

E = (constant) x exp (— 2rjt). 

It follows that: 

log-E = — 2r/t + log (constant), 
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and taking derivatives: 

dE/dt = -2rjE. (2.12) 
If we combine Equations 12.101 12.111 and 12.121 we have: 

dE/dt = E oc -T{v)I. (2.13) 

We may attempt to measure solar-cycle variations in A, H, P, E and E. The pa- 
rameters that are usually extracted directly by the observers are the widths (A) and 
heights (H) of the peaks in the frequency power spectrum. The mode powers, P, are 
then estimated from 

P=^AH (2.14) 

(i.e., area under Lorentzian). The E may be computed from the P using model computed 
inertias, /. We note that fractional changes in E will be the same as those in P (assuming 
the I do not show any significant changes). Finally, E follows from 

E = -irEA = -ttIPA. (2.15) 

If solar-cycle variations are uncovered in the above parameters, what do they imply for 
the underlying cause (or causes)? If we think in terms of the basic oscillator parameters - 
the damping rate r] and the forcing function J-(v) - Equations l2.81 [231 12. 10112. Ill and l2. 131 
then imply the following: 

SAv/Av = 6r]/r), (2.16) 
SH/H = 6T(u)/T(u) -26r)/r), (2.17) 
5P/P = SE/E = 5F{v)/F{v) - 5r)/r), (2.18) 

and 

SE/E = ST(v)/T(v). (2.19) 

Observations of p modes at low I are consistent with SE/E — 0, which in turn implies 
that 8J-(v) / J-(v) = 0: the acoustic forcing of the low-/ modes remains, on average, 
constant over the cycle. Any changes that are observed in the other parameters must 
therefore arise from Sr\/r\ alone, and we have that: 

SAv/Av = -SP/P = -SH/2H, (2.20) 

The observed changes in A, P and H follow these ratios. We may therefore infer that 
changes to these parameters are in all likelihood the result of changes to the damping 
of the low-/ modes. Houdek et al. (2001) noted that the high Rayleigh number of the 
magnetic field could modify the preferred horizontal length-scale of the convection, which 
in turn could modify the damping rates. Their numerical calculations of the implied 
changes for low-/ modes show reasonable agreement with the observations. 



3. Subtle signatures of stellar activity and stellar-cycle variations 

Asteroseismic observations of stars will, for the foreseeable future, be limited to ob- 
taining data on low-/ modes because of the cancellation effects resulting from unresolved 
observations of stellar discs. Moreover, the angle of inclination, i, offered by a star also 
determines which of the low-/ components will have non-negligible visibility in the data. 
The visibility of the radial (/ = 0) modes is not affected by i, but for the non-radial 
modes the visibility, £{i)i m (in power) of a given (/, m) component is assumed to go like 
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(e.g., Gizon & Solanki 2003): 

where P ; ' m are Legendre polynomials. The fact that some components may in effect 
be missing from the data has important implications for the analysis. There are some 
obvious implications, e.g., it will not be possible to measure frequency splittings of non- 
radial modes if i is close to 90 degrees - so that the star is viewed rotation-pole on 
- since only the zonal components (to = 0) will have non-negligible visibility. In this 
section we highlight a more subtle aspect, which relates to the fact that estimates of the 
frequencies of non-radial modes are affected by the interplay between the effects of near- 
surface magnetic activity and the relative visibility (or absence) of the mode components. 
This has consequences for inferences made on solar-type stars. We first introduce and 
explain the problems using "Sun-as-a-star" helioseismic observations as a test case, before 
considering the wider implications for asteroseismology. We return later, in Section [4] to 
discuss the implications for estimation of average frequency shifts of solar-type stars. 

3.1. Explanation and description of frequency bias 

Extant Sun-as-a-star observations are made from a perspective in which the plane of the 
rotation axis of the Sun is nearly perpendicular to the line-of-sight direction (i.e., i is 
always close to 90 degrees) . This means that only components with / + m even have non- 
negligible visibility. As we shall explain below, the impact of surface magnetic activity 
on the arrangement of the components in frequency means it may not then be possible 
to estimate the true frequency centroid of the multiplet. The frequency centroids carry 
information on the spherically symmetric component of the internal structure, and are 
the input data that are required, for example, for hydrostatic structure inversions. 

In the complete absence of the near-surface activity, the / + to odd components which 
are "missing" from the Sun-a-as-a-star data would be an irrelevance. All mode compo- 
nents would be arranged symmetrically in frequency, meaning centroids could be esti- 
mated accurately from the subset of visible components. A near-symmetric arrangement 
is found at the epochs of modern solar-cycle minima (Chaplin et al. 2004b). However, 
when the observations span a period having medium to high levels of activity - as a 
long dataset by necessity must - the arrangement of components will no longer be sym- 
metric. The frequencies given by fitting the Sun-as-a-star data will then differ from the 
true centroids by an amount that is sensitive to I. The I dependence arises because in 
the Sun-as-a-star data modes of different I comprise visible components having different 
combinations of I and to; and these different combinations show different responses (in 
amplitude and phase) to the spatially non- homogeneous surface activity (i.e., as observed 
in the active bands of latitude for solar- type stars) . 

These effects are illustrated in Figs. |6] and [7] Fig. [6] is a schematic plot of the ar- 
rangement in frequency of the components of an I = 2 mode, under circumstances where 
activity levels are assumed to be negligible. The top panel shows the arrangement in 
frequency of the various to (as labelled on the ordinate) , assuming a rotational frequency 
splitting of 0.4 /xHz between adjacent to. The dotted line marks the frequency centroid. 
The middle panel shows schematically both the placement in frequency (abscissa) and 
relative visibility in the frequency power spectrum (height on ordinate) of the three mode 
components that would have non-negligible visibility in the Sun-as-a-star data. The bot- 
tom panel shows the Lorentzian peak profiles expected of a mode at the centre of the 
p-mode spectrum (where peak linewidths are around 1 /iHz). 

Fig. [7] shows the situation at an epoch coinciding with high levels of solar activity. 
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Figure 6. Top panel: Placement in frequency, relative to the frequency centroid of all 21 + 1 
components (abscissa) , of the various m (labelled on ordinate) of an / = 2 mode during an epoch 
of negligible solar surface activity. The solid vertical line marks the location of the centroid. 
Middle panel: placement in frequency (abscissa) and relative visibility in the frequency power 
spectrum (height on ordinate) of the three mode components that would have non-negligible 
visibility in the Sun-as-a-star data. Bottom panel: Lorentzian peak profiles expected of an I = 2 
mode at the centre of the p-mode spectrum (where peak linewidths are around 1 fiRz). 

This activity is most prominent in the active latitudes. The |m| = 2 components are 
more sensitive to acoustic perturbations due to activity that arise in these latitudes 
than is the m = component (as a result of the spatial distribution of the relevant 
spherical harmonic components) . The result is a distortion of the symmetric arrangement 
in frequency of the multiplet, because the different |m| experience a frequency shift of 
different size. 

The non-symmetric arrangement of the components is very evident in Fig. [7j Two 
lines have been plotted in the middle panel, which were not shown in Fig. |6l The solid 
line marks the (unweighted) average frequency of the components showing non-negligible 
visibility. It is clearly offset, or biased, by about 0.07 //Hz, from the centroid (dotted line). 
The dashed line marks the actual, expected location of the "Sun-as-a-star" frequency 
(offset from the centroid by more than 0.1 /iHz). It differs from the unweighted average 
of the three visible components because the "peak bagging" procedures used to estimate 
the mode frequencies are influenced by the relative heights of the different m. Here, the 
outer \m\ = 2 components carry greater weight in determining the fitted I = 2 frequency 
than does the weaker m = component, thereby shifting the estimated frequency from 
the unweighted average. 

It can be quite hard to estimate, to any reasonable accuracy and precision, the loca- 
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Figure 7. Same as Fig. [6] but during an epoch with activity levels corresponding to the max- 
imum of a typical (modern) solar cycle. The solid and dashed lines in the middle panel mark 
estimates of the I = 2 frequency that might be derived from the Sun-as-a-star data (see text). 



tions in frequency of each component individually (which, as noted above, would help to 
reduce the offset between the estimated frequency and the true centroid because it would 
be possible to compute the unweighted average frequency). Frequency estimation must 
contend with the fact that within the non-radial multiplets the various m lie in such 
close frequency proximity to one another that suitable models which seek to describe the 
characteristics of the m present must usually be fitted to the components simultaneously. 

Is there any way to get around the problem, and in some way "correct" the Sun-as- 
a-star frequencies for the frequency bias? Two methods have been applied successfully. 
The first relies on the availability of contemporaneous resolved-Sun helioseismic data 
- which are sensitive to all m components - to in effect map the strength and spatial 
distribution of the acoustic asphericity arising from the surface activity over the epoch 
in question. Chaplin et al. (2004c) and Appourchaux & Chaplin (2007) show how to 
make the correction, using the so-called "even a coefficients" from fits to the resolved- 
Sun frequencies (the a coefficients are discussed later in this section). This approach is 
of course currently not viable for other solar-type stars. Fortunately, the second method 
in principle is: here, the same dataset that is used to produce the frequency estimates 
is also used to help calibrate the corrections that are needed to "clean" the frequencies. 
Implicit in the procedure is the assumption that at epochs coinciding with cycle minima, 
any frequency asymmetry in the mode multiplets is negligible. This assumption turns 
out to be reasonable for the Sun, as noted previously (Chaplin et al. 2004b). 

The second method is described at length in Broomhall et al. (2009). In sum, the long 
dataset is divided into smaller subsets, and the frequencies of each subset are estimated. 
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Tried and trusted procedures are then used to measure the frequency shifts of different 
modes from one subset to another. The frequency shift of any given mode may be mea- 
sured relative to its average frequency over the long dataset, or, for example, its frequency 
in the subsets (or subset) which coincide(s) with the cycle minimum (the choice has sub- 
tle consequences for propagation of the frequency uncertainties, e.g., see Chaplin et al. 
2007b). It is common practice, at least for the low-^ data, to average shifts over several 
modes to reduce uncertainties on the shifts. By making an average over many orders in 
n, one can track the average frequency shift from subset to subset. By making an average 
over a few (say two, or three) orders in n it is possible to constrain the dependence of 
the shifts on frequency (or n). 

In the next stage of the procedure the frequency shifts are parametrized. Sun-as-a-star 
procedures have parametrized the shifts as a linear function of one of the global proxies 
of magnetic activity, the 10.7-cm radio flux (Tapping & De Tracey 1990) being a favoured 
choice. The mode-averaged frequency shifts are fitted to a linear function of the 10.7- 
cm flux (averaged over the same periods as the seismic data). This linear model then 
suffices to predict the expected average frequency shift for any epoch, given the 10.7- 
cm flux for that period. Expected shifts for individual modes are then computed using 
a best-fitting polynomial, which describes the measured frequency shifts as a function 
of frequency, relative to the average frequency shift. The predicted frequency shifts may 
then be subtracted from the raw frequencies (the latter measured from the full timeseries) 
to yield a set of corrected frequencies, commensurate with negligible surface activity. The 
procedure yields frequencies that are not subject to the bias discussed above, since that 
bias is assumed, by definition, to be missing at low activity. 

Application of the procedure, verbatim, to other stars would of course demand that 
contemporaneous activity data be available (e.g., chromospheric H & K data). Such data 
may be hard to obtain in sufficient quantities on a large number of stars. One might 
instead use measures of variability in the timeseries - arising from rotational modulation 
of surface activity - to form a suitable proxy (Basri et al. 2010). It may instead be 
possible to parameterize the frequency shifts by fitting non- linear functions in time (i.e., 
to describe the observed shifts as functions of time), hence circumventing the need to 
have complementary activity data. 

How large might the frequency bias be for other stars? The magnitude of the bias is 
governed to first order by the size of the acoustic asphericity, arising from the magnetic 
activity, and the angle of inclination of the star. Second-order effects come in through the 
interplay between the mode properties and the peak-bagging procedures used to estimate 
the frequencies (e.g., see discussion in Ballot et al. 2008 on correlations between different 
parameters) . 

Fig. [8] shows the predicted bias from the first-order effects. Here, we assumed an acous- 
tic asphericity similar to that observed on the Sun, as experienced by I = 1 and I = 2 
modes at the centre of the solar p-mode spectrum. The predictions therefore in principle 
show the bias that would be expected in the most prominent low-^ modes if the Sun were 
to be observed at different angles of inclination, i. Two predictions are plotted in each 
panel. We first explain how we derived the solid-line prediction (using the I — 2 case as 
an example). 

The full set of m frequencies of a mode may be described by a polynomial expansion 
of the form 



j=2l 




(3.2) 
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where v ccn is the frequency centroid of the mode, and V\ (m) are polynomials related 
to the Clebsch-Gordan coefficients (Ritzwoller & Laveley 1991). An expansion of this 
type is commonly used to fit the frequencies observed in resolved-Sun data, where all m 
components are available. The even a,j describe perturbations that are non-spherically 
symmetric in nature, e.g., the acoustic asphericity from activity, while the odd cij coeffi- 
cients describe spherically symmetric contributions to the frequency splittings (rotation) . 
Here, for simplicity, we assume that for the low-? modes the a 2 term dominates other 
terms in the description of the asphericity, and that the a\ term dominates other terms 
in the description of the rotation. The required V\ (to) are: 

Tl(m)=m/l (3.3) 

and 

a 6m 2 -2/(Z + l) 

Vl{m)= 6P-M0 + 1) - (3 - 4) 

Wc may then write the frequencies of each of the components explicitly as: 

V n 20 = V CC n ~ 2o 2 , (3.5) 

and 

V n 1\ = Vccn + a\ — 0,2, 

v n i _i = v csa — a\ — 02, (3.6) 

and 

v n 22 = Vccn + 2ai + 2a 2 , 

v n 2 -2 = ^cn - 2ai + 2a 2 . (3.7) 

The solid line prediction in Fig. [8] assumes that the estimated frequency of the mode 
is the weighted combination of the individual frequencies, where the weights are propor- 
tional to the relative heights of the components, £(i)i m (see Equation 13. ip . We weight 
in this way because the estimated frequency is influenced by the relative heights of the 
different m. We then have: 



/m=+2 \ /m=+2 

Ai2 = V 



£(i)2m X V n2m / £{i)2m ■ (3-8) 



m— — 2 / \m— — 2 



The solid lines in Fig. [8] show v n i ~ v ccn , i.e., the predicted frequency bias, as a function 
of i. To make the prediction we adopt a 2 = 0.1 /iHz (typical value at solar maximum). 
(The results do not depend upon ai, which cancels when v n i is computed.) 

Results on Sun-as-a-star data, and simulated asteroseismic data on solar-type stars, 
suggest that the weights are actually non-linear functions in £ (i)i m (e.g., see Chaplin ct 
al. 2004b). A better description is one for which the weights are proportional to £(i)i m 
raised to some positive power (four is probably a reasonable value). Here, we therefore 
also make predictions based upon 



Vn2 



>< "n2m / £ £«2rn , (3-9) 



which gives the dashed lines in Fig. HI 
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Figure 8. Expected frequency bias for Sun-as-a-star observations made for different angles of 
inclination, i (see text for explanations of solid and dashed lines). 

Equations 13.81 and 13.91 may be written in more general form (i.e., for any given I) as: 

(m=+l \ /m=+l \ 

E £«7 m x*W / £ £ Wm > (3-10) 
m— — / / xm=—l / 

where 7 is a constant. 

The results indicate that the bias is least severe at intermediate values of i. The worst 
case is when i = degrees. The non-linear weighting in £{i)i m also changes slightly the 
functional form of the bias in i, e.g., in the case of I = 1 it extends the range of i over 
which the bias is more important. 

When we obtain long datasets on solar-type stars, frequency uncertainties will drop 
to levels well below 0.1 /iHz. In this context, the bias at some angles is by no means 
insignificant. Moreover, the bias would of course be larger on a star with higher-than- 
solar acoustic asphericity, i.e., one might quite reasonably expect to encounter solar- type 
stars where the bias may be a significant fraction of 1 /iHz or more. In sum, this source 
of bias in the frequencies will need to be taken into account in the not-too-distant future, 
if we are to avoid errors on the inferences drawn on target stars. 

3.2. Impact on frequency separation ratios 

Roxburgh & Vorontsov (2003) proposed the use for asteroseismology of the ratio of 
separations in frequency between \ow-l p-modes: Frequency separation ratios are formed 
from the large and small frequency separations of the p modes, 

Afi(n) = v n i - Vn-1,1, (3-11) 



and 

du+i{ri) = v n ,i - v„-ij+2- (3-12) 
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Data on modes in the range < I < 3 may be used to make ratios, according to: 

d 02 (n) d 13 (n) 

rM = A^ry ri3(n) = AMn+iy (3 - 13) 

Roxburgh & Vorontsov used an asymptotic method to show that the ratios are very in- 
sensitive to conditions in the near-surface layers of a solar-type star. This is an attractive 
property since the outer layers are uncertain and hard to model accurately. As such, the 
ratios in principle offer a clean diagnostic of the internal properties of a star (see also Otf 
Floranes et al. 2005 and Roxburgh 2005). They are nevertheless potentially sensitive to 
the acoustic asphericity - as was pointed out by Oti Floranes et al. (2005), and discussed 
at length by Chaplin et al. (2005) - which is clearly undesirable if one wishes to suppress 
any signatures from the near-surface layers. 

That the ratios are sensitive to the solar cycle may be understood as follows. For the 
Sun-as-a-star data, measured frequency shifts of the I = 2 modes are larger than those of 
the neighbouring / = modes (because the frequencies of these modes are dominated by 
the \m\ = 2 components, which are more sensitive to the acoustic perturbations arising 
in the active latitudes than are the I = modes). The measured di 1+2(11) will therefore 
carry a residual signature of the solar cycle, being smaller at higher levels of activity. 
Residual solar-cycle changes in Avi(n) are, fractionally, much smaller (but nevertheless 
measurable; see Broomhall et al. 2011), so that in the case of, for example, r^n): 

5r 2(n) /r 2(n) ~ 5d 2(n)/do2(n). (3-14) 

For the Sun-as-a-star data, 6do2(n) / do2(ri) ~ —0.01 for modes at the centre of the p- 
mode spectrum observed at solar maximum, hence, we would expect the ratios to decrease 
between solar minimum and solar maximum by about this amount. Fig.|9]plots measured 
fractional differences in ro2, as averaged over the most prominent low-Z modes. The 
differences were computed between two BiSON Sun-as-a-star datasets, one formed of 
from around solar maximum, the other formed of data from around solar minimum. 

The potential impact on observations of other solar-type stars will depend upon both 
the nature of the near-surface magnetic fields, and also the angles of inclination offered 
by the stars. In Section 0] below, we shall develop a simple model of the frequency shifts 
that will allow us to make these predictions for other stars. We therefore return in that 
section to the frequency separation ratios. 



3.3. Impact on shapes of resonant 

What effect will stellar-cycle shifts in frequency through a long timeseries have on the 
underlying shapes of the p-mode peaks, when those peaks are observed in frequency- 
power spectra made from the full timeseries? Given a significant shift, one would expect 
the Lorentzian-like peak shapes to be distorted, with the mode power being spread out in 
frequency thereby flattening the profiles. Whether or not a frequency shift 5v is significant 
in this context is determined by the ratio 

e = 8v/A, (3.15) 

i.e., the ratio of the shift to the linewidth of the mode. The larger this ratio, the more 
distorted a mode peak will be. 

Chaplin et al. (2008b) derived analytical descriptions of the distorted peak profiles. 
Since for solar-type stars we expect the modes to be excited and damped on timescales 
much shorter than that on which any significant change of the mode frequencies is ob- 
served, the profiles may be assumed to correspond to the average of all the instantaneous 
profiles taken at any time t within the full period of observation T. Each instantaneous 
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Figure 9. Fractional change in ro2 - as averaged over the most prominent low-i modes - between 
two BiSON Sun-as-a-star datasets, one around solar maximum and one around solar minimum 
(differences in sense maximum minus minimum). (From Chaplin et al. 2005.) 



profile may be described as, for example, a Lorentzian with a central frequency v(t)t, 
such that the time-averaged profile is 



{P(y)) = - 



H 



dt, 



1 



v-v(t) 
A/2 



(3.16) 



where the angled brackets indicate an average over time, and H and A are the mode 
height (maximum power spectral density) and linewidth, respectively. 

Chaplin et al. (2008b) derived profiles given by two functions describing the frequency 
shifts in time: first, the simplest possible function, this being a linear variation over 
time (which would be relevant for describing the profiles given by observations made 
on the steepest parts of the rising or falling phases of a stellar cycle); and second, a 
co-sinusoidal variation, to mimic a full stellar cycle. They neglected the effects of the 
solar-cycle variations in mode power, height and width. 

A simple linear variation in time may be described by 

v(t) = vo + 8v^ (3.17) 

where v$ is the unperturbed frequency, and the frequency is shifted by an amount 8v 
from the start (t — 0) to the end (t — T) of the timeseries. Substitution of Equation 13. 171 
into Equation [3461 followed by solution of the integral, gives the predicted mode profile: 



2c 



+ X 2 

where e is the shift-to-linewidth ratio (Equation I3.15|) . and 

v _ v - (i/„ 4- 5v/2) 
X - Xf2 ■ 



(3.18) 



(3.19) 
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The top panel of Figure \TU\ shows profiles given by Equation 13.181 The unperturbed 
profile (solid line) is for a mode having an unperturbed frequency of vq = 3000 /iHz, 
an unperturbed linewidth of A = 1 //Hz, and an unperturbed height of H = 100 units. 
The other curves show the profiles that result when the frequency shift, 6v, is: 0.15 //Hz 
(dotted line); 0.40 /iHz (dashed line); 1.50 /iHz (dot-dashed line); and 3.0 ^Hz (dot-dot- 
dot-dashed line). Since A = 1 /iHz, the 8v also correspond to the shift-to-linewidth 
ratios, e. To put the values in context, we once more recall that low-^ solar p modes at 

3000 /iHz, which also have width rs 1 /iHz, and show a frequency shift of approximately 
0.40 /iHz from the minimum to the maximum of the solar activity cycle. 

Only at the two largest shifts (dot-dashed and dot-dot-dot-dashed lines) do the profiles 
depart appreciably from the Lorentzian form. Closer inspection of the profiles does reveal 
some modest distortion at the two, smaller, Sun-like shifts. These have e = 0.15 and 0.40 
respectively. 

A co-sinusoidal variation of the mode frequency may be described by 



u{t) = v + — 



1 — cos 



2vrf 
~p 

1 eye 



(3.20) 



where vo is again the unperturbed frequency (i.e., the frequency at minimum activity), 
while 8v is now the full amplitude (from minimum to maximum) of the cyclic frequency 
shift (not the total shift, as in the linear model). When the length of observation, T, 
equals the cycle period P cyc (or one-half of the period), the average profile resulting from 
Equation l3~2"0l is: 

( p o)> = tHwtt (3 - 21) 

y/1 - e 2 .F{L 1 ,L 2 ) 



where 

F{L U L 2 ) 

Li 
L 



4LiL 2 



H 



l + {X~ef : 



H 



2 l + (X + ef' 
iz-fo + f) 
X - Xf2 ' 

Since the frequency spends more time around its maximum and minimum values, power 
near these extreme frequencies will have more weight in the time-averaged profile, giving 
the profile a double-humped appearance. This is reflected in the analytical expression 
through the two Lorentzians L\ and L 2 . It is obvious from Equation 13.211 that when 
f« 1 the profile tends to a single Lorentzian. 

The bottom panel of Figure [T0l shows predicted profiles from Equation l3.21ll|i . assuming 



observations made over a complete activity cycle, and with the same shifts 8v as were 
applied in the linear-model case (top panel of Figure [TU)) . As for the simpler linear 
variation, it is only at the two largest 8v that the profiles depart appreciably from the 
unperturbed (Lorentzian) form, here showing the predicted "humps" at the extreme 
frequencies of the cycle. However, closer inspection again reveals some distortion of the 

f Versions of Equations 13. 181 and 13. 21 1 which allow for the small amounts of peak asymmetry 
observed in real solar p modes, are presented in Chaplin et al. (2008). We note here that this 
asymmetry is negligible for modes at the centre of the solar p-mode spectrum. 
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Figure 10. Top panel: Peak profiles expected for a single mode of width A = 1 piHz in the 
frequency-power spectrum of a time series within which the frequency was varied in a linear 
manner by total amount 8v. Various linestyles are: no shift (solid line); &v = 0.15 fiRz (dotted 
line); 0.40 /iHz (dashed line); 1.50 fiHz (dot-dashed line); and 3.0 /iHz (dot-dot-dot-dashed line). 
Bottom panel: Peak profiles expected for a single mode of the same width, but where the 
frequency has instead been subjected to a co-sinusoidal variation in time. The timeseries is 
assumed to have length equal to one cycle period; while the amplitude of the cycle is 5v/2 
(giving a total minimum-to- maximum shift in frequency of 5v). Linestyles are as per the upper 
panel. (From Chaplin et al. 2008b.) 

Lorentzian shapes at the small Sun-like shifts. For a given shift, this distortion appears 
to be slightly larger than in the simpler, linear case. 

While the distortion seen at the Sun- like shifts is evidently modest, Chaplin et al. 
(2008b) showed that it is nevertheless just detectable in long, multi-year Sun-as-a-star 
datasets. By fitting modes to the usual Lorentzian-like models - which do not allow for 
the distortion - rather than modified models, like the ones shown above, Chaplin et al. 
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showed that an overestimation (underestimation) of the linewidth (height) parameter 
results. This bias is estimated to be of size comparable to the observational uncertainties 
given by datasets of length several years. Bias in the frequency parameter is much less 
of an issue. 

The distortion may of course be more important for asteroseismic datasets on some 
solar- type stars, e.g., those for which the ratio of the stellar-cycle frequency shifts to the 
mode linewidths is larger than for the Sun. The shifts need only be about twice as strong 
as those on the Sun before significant distortion of the peaks results. Visible distortion 
of the mode profiles in asteroseismic data may as such provide an initial diagnostic of 
strong stellar-cycle signatures over the duration of the observations. 

4. Estimated frequency shifts, and inference on spatial distribution of 
surface activity 

We have seen that near-surface activity that is distributed in a non-spherically sym- 
metric manner (i.e., the aforementioned acoustic asphericity) will give rise to differences 
in the magnitudes of the frequency shifts of modes of different (I, \m\). For a scenario like 
that on the Sun, where bands of strong magnetic activity are located at lower latitudes, 
it will be the sectoral components of the non-radial modes that show the largest shifts. 
As was noted in Section 13.11 these components, for which I = |m|, have more sensitivity 
to changes at lower latitudes than do their zonal counterparts. It is possible to measure 
differences in the sizes of the frequency shifts of the different low-/ modes (e.g., see Chap- 
lin et al. 2004a; Jimenez-Reyes 2004). If such measurements could be made on other 
stars, might it not then be possible to make some inference on the surface distribution 
in latitude of the activity, from measurement of the relative sizes of the frequency shifts 
(assuming, as mentioned elsewhere, that it is near-surface perturbations that are the 
dominant cause of the observed shifts)? 

This question was considered by Chaplin et al. (2007a). They used a very simple model 
of the spatial distribution of the surface activity on a solar-type star, in which the activity 
was given a uniform amplitude between some lower and upper bounds in latitude, A. As 
should be evident from what follows, it would be very straightforward to incorporate 
a more sophisticated model of the activity. They used this model to predict expected 
frequency shifts of modes of different (I, \m\). Here, we reproduce their model, and use 
it to show explicitly how it might be used with results on observed frequency shifts to 
make inference on the active latitudes of a star. Before that, we use its predictions to 
show how estimated average frequency shifts, and estimated frequency separation ratios, 
of solar-type stars are affected by the angle of inclination. 

Implicit in the calculation of the expected frequency shifts is the assumption that 
contributions to the fractional change in the sound speed are non-zero only very close 
to the surface. Consideration of the radial dependence of the sound speed is neglected 
(again, to simplify the description). The expected shifts, 5vi m are then proportional to 
(e.g., Moreno-Insertis & Solanki 2000) 



where the P/™(cos#) are Legendre polynomials, and B(9) describes the distribution in 
co-latitude of the activity. (Note that A = 90 degrees — 6. ) Chaplin et al. calibrated the 
model so that for a Sun-like configuration the predicted Svoo was equal to the observed 
average shift of the solar I = modes. This calibration was subsumed within the B(8), 
as discussed below. For all surface configurations tested, the activity at cycle minimum 




(4.1) 
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Figure 1 1 . Frequency shifts given by a simple model of stellar surface magnetic activity (mag- 
netic flux). Linestyles, which are the same in all panels, show results for different (Z,|m|) com- 
ponents as follows: solid for (0,0); dotted for (1,0); short dashes for (1,1); dot-dashed for (2,0); 
dot-dot-dot-dashed for (2,1); and long dashes for (2,2). Top panels: shifts calculated from models 
based on Equation 14.21 (left-hand panel) and Equation 14.31 (right-hand panel). Middle panels: 
residuals given by subtracting the (0,0) shifts from the other component shifts. Bottom panel: 
ratio of component shifts to (0,0) shifts (these ratios are the same for both sequences of models). 
(From Chaplin et al. 2007a.) 
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was set to B(9) = 0. Values of 8vi m determined at the simulated cycle maxima therefore 
corresponded to the sought-for stellar cycle shifts. 

In a first sequence of models B{9) was set to some constant value at cycle maximum. 
This value was the same for all model configurations. The only parameter that was 
changed from model to model, and therefore the only factor that could alter the observed 
shifts, was the maximum latitude of the activity, A max . The minimum latitude was fixed 
for all models at A m j n = 5 degrees. In summary: 

const for A min < |A| < A max , 
BU-h { (4.2) 
otherwise. 

In the scenario above, the total surface magnetic flux is proportional to A max — A m i n , 
and therefore varied by over an order of magnitude for the range of computations made 
(10 < A max < 70 degrees). Chaplin et al. also made a second sequence of models in which 
the total flux was conserved for all values of A max . This model is described by: 

const/ (A max - A min ) for A min < |A| < A max , 
BW) = { (4.3) 
otherwise. 

Chaplin et al. used results from Knaack et al. (2001) to help guide their calibration 
of Equations 14.21 and 14.31 Knaack et al. noted that observations of the surface activity 
during the cycle 22 maximum showed sunspots confined to bands from 5 to 30 degrees, 
and faculae in bands from 5 to 40 degrees. On the Sun, faculae occupy a much larger 
surface area than the spots (see also De Toma et al. 2004). The constant in Equations l4.2l 
and 14.31 was therefore calibrated so that when A max = 40 degrees, the calculated shift of 
the radial modes was 8vm ~ 0.3 /iHz. This is the value observed for the most prominent 
low-Z radial modes of the Sun. 

The top two panels of Fig. [TT1 show the calibrated shifts from both sequences of models, 
based on Equation 14.21 (left-hand panel) and Equation 14.31 (right-hand panel). The mid- 
dle panels show residuals given by subtracting the (I, |m|)=(0,0) shifts from the other 
component shifts. The bottom panel shows the ratio of the component shifts and the 
(0,0) shifts: since B(9) = const, for both sequences of models, the ratios are also the 
same (hence only one plot). 

Let us first consider the impact of the angle of inclination, i, on estimated average 
frequency shifts for solar-type stars. We again adopt the approach of Section 13.11 to 
take into account the effect of the relative visibility of components within any given / 
(Equations 13. 81 to 13. 10p . We take the frequency shifts 8vi m predicted by the simple model 
in this section, for different angles of inclination, i (as computed assuming the solar 
calibration). Since the 8vi m are functions of i, we write Sui m (i). We weight the shifts of 
each component within a multiplet, to give the effective (weighted average) shift at that 

(m=+l \ /m=+l \ 

£ £(i)7 ro x^J/( ( 4 - 4 ) 
m— — l / \jn— — l / 

Since one may then choose to average results over different I, to reduce errors, we compute 
a final, weighted average shift according to: 

(^)=^E^ X ^-J/(E^)' (4.5) 
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Figure 12. Predicted frequency shift, averaged over I = 0, 1 and 2 modes, as a function of 
angle of inclination, i, assuming a solar-like spatial distribution of surface activity (see text for 
details) . 



using the predictions for / = 0, 1 and 2. Here, the £[ are the relative total / visibilities, 
which we take to be 1, 1.5 and 0.5 (for I = 0, 1 and 2, respectively). 

Fig- HUpkots the (8v) as a function of i, for 7 = 1 (solid line) and 7 = 4 (dashed lines). 
When i is low, the m = components have the highest visibility in the I = 1 and I — 2 
modes, which reduces the size of the average frequency shift (since these components 
show smaller shifts than the other components when A max = 40 degrees) . The plot shows 
clearly that for a solar-like spatial distribution of surface activity, estimated average 
frequency shifts could differ by up to a factor of two, depending upon the angle i. It 
will clearly be important to have good constraints on i in order to properly interpret 
measured frequency shifts of solar-type stars (e.g., see Ballot et al. 2006, 2008). 

We also use predictions from Equation l4.4l above to show the impact of i on the observed 
frequency separation ratios. Fig. [13] plots the implied fractional change, between solar 
minimum and maximum, in the frequency separation ratios ro2 of the most prominent 
solar p modes. Recall from Section 13.21 that results from BiSON data give a fractional 
change of —0.01; this is consistent with the prediction shown here for i = 90 degrees. 
Fig. [13] shows that for the Sun the largest offset occurs at i = degrees, when the 
predicted fractional change is 0.02 (positive). Since we should expect to observe solar- 
type stars with larger acoustic asphericity than the Sun, it is possible that bias in the 
frequency separation ratios could reach levels of several per cent. 

Let us now use the frequency shift ratios plotted in the bottom panel of Fig. 111! 
together with results on measured frequency shifts of low-Z solar p modes, to make an 
estimate of A max for the Sun. By using the ratios as opposed to the absolute values of the 
shifts we remove any dependence of the result on the absolute calibration of the model 
shifts. Let us denote the shift ratios by fc; m = 8vi m /8vQQ. From the BiSON Sun-as-a-star 
data in Chaplin et al. (2004a), we may infer observed ratios of fen = 1.23 ± 0.14 and 
k-22 = 1.45±016. (The observed I = 2 Sun-as-a-star shifts will include a contribution from 
the (2,0) components, however the contribution is relatively weak, and the observed shifts 




Figure 13. Predicted fractional change (from solar minimum to solar maximum) in frequency 
separation ratios ro2 of the most prominent solar p modes, as a function of angle of inclination, 
i, assuming a solar-like spatial distribution of surface activity (see text for details). 



may be regarded, to very good approximation, as being those of the (2,2) components.) 
To infer A max , we use the appropriate curves from Fig. 111! Uncertainties in A max follow 
from 

<5A max = ^ max Jfc im , (4.6) 

where 8k i m are the observed uncertainties in the fc; m , and AA max /Afc/ m are the gradients 
of the respective curves in the vicinity of ki m . From the I = 1 shift ratio we find A max = 
43 ± 16 degrees, and from the I = 2 shift ratio we find A max = 34 ± 10 degrees. These 
values are of course reasonable estimates for the Sun. 

One may adopt a similar procedure using asteroseismic data on other solar-type stars. 
One would again need to have good constraints on i, which may be obtained from peak- 
bagging of the modes in the frequency-power spectrum. 



5. How long will it take to detect evolutionary changes? 

We finish by thinking longer-term, and ask whether it might be possible to measure 
evolutionary changes in solar p-mode frequencies, using low-^ data (the choice of data 
again made with other solar- type stars in mind). 

First, consider the precision in the estimated frequencies. We take the estimated un- 
certainties, (T n i, from Broomhall et al. (2009), who measured low-Z frequencies in 8640 d 
(23.7 yr) of BiSON Sun-as-a-star data. The published list included 81 frequencies, cov- 
ering I — to 3. An estimate of the combined precision in all the measured frequencies, 
(a), is given by: 

( <T )=fe 1 /^ ^l.lnHz. (5.1) 
We may then estimate the combined precision, a At, for any dataset length, At (in yr) 
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Figure 14. Combined frequency precision, a At, for low-Z observations, extrapolated to 
different dataset lengths, At (see text for more details). 



from: 



o~At 



( 8640 \ 1/2 
V~365~ J 



a) At- 1/2 ~ 5.4At- x / a nHiz, 



(5.2) 



since the frequency uncertainties of the observed modes are expected to scale with the 
square root of time. Fig. [TJ] plots a At as a function of dataset length, At. 

How do the a At compare to the expected frequency shifts from evolutionary effects? 
Those shifts should, to good approximation, be dominated by the slow expansion in 
radius of the Sun as it evolves on the main sequence. Since, to first order, the frequencies 
scale like the square root of the mean density of the star, we have: 



3 /SR 
2 l~R 



This implies that 











2R 


\St) 



(5.3) 



(5.4) 



Fig.[15]plots l/v(8v/v), as determined from the frequencies of several stellar evolutionary 
models (Yale-Yonsei models) made with solar mass and solar composition, but with 
varying ages ranging ±0.02 Gyr about the accepted solar age. The different linestyles 
show results for different / (/ = as a thin solid line, I — 1 as a dotted line, I = 2 as a 
dashed line, and I — 3 as a dot-dashed line). The thick grey line plots — (3/2R)(5R/ R), as 
determined from the radii of the stellar models. It has a value —5.3 x 10 -11 yr -1 , which 
agrees quite well with the estimated gradient from the model-computed frequencies, 
verifying that the frequencies scale in an approximately homologous fashion. 

We may then estimate to good approximation the expected evolutionary frequency 
shift Av n i of a mode over a length of time At (in yr): 



Av„ 



-5.3 x lQ- n i/ nl At 



(5.5) 
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Figure 15. v{Sv/v), as determined from the frequencies of several stellar evolutionary models 
made with solar mass and solar composition, but with varying ages, ranging ±0.02 Gyr about 
the accepted solar age. Linestyles show data for I — (thin solid line), I — 1 (dotted line), 
I = 2 (dashed line) and I = 3 (dot-dashed line). The thick grey line plots —(3/2R)(5R/R), as 
determined from the radii of the stellar models. 

For a mode at the centre of the solar p-mode spectrum, having a frequency w 3000 /iHz, 
we obtain —1.6 x 10~ 4 nHz in 1 yr, and —1.6 x 10~ 2 nHz in 100 yr. These expected shifts 
are significantly smaller than the combined uncertainties expected for each length, being 
just over 5 nHz in lyr, and about 0.5 nHz in 100 yr. The shift and combined uncertainty 
have similar sizes only when At w 1100 yr. 

The author would like to express his thanks to P. Palle and colleagues at the IAC, for 
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